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. Introduction : (a) observation 

Fibrous Thermal Protection Systems (TPS); e.g. Stardust and PICA 
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First step for the 
porous medium 
approach : oxidation 

[1 ] M. Stackpoole et al., Post-Flight Evaluation of Stardust Sample 
Return Capsule Forebody Heatshield Material, AIAA 2008-1202 


SEM 3 and 4 suggest the 
occurrence of some volume 
phenomena in the char layer: 

- Oxidation (oxygen from the 
atmosphere) 

- Sublimation 

- Mechanical erosion of the 
matrix 



SEM micrographs h) 
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. Introduction : (b) comparison of 2 materials 

Equilibrium chemistry vs. Finite-rate chemistry 


Same fibrous preform, chemical composition, overall density 
A : dense matrix layer around the fibers B : Expanded, low density pore-filling matrix 



“Surface ablation” model (as described by Kendall et al., NASA CR-1060, 1968) 
Equilibrium chemistry -> Only the chemical composition is important 

- model for A = model for B (in a control volume above the effective surface) 
“Ablation-zone model” model 

- Finite-rate chemistry -> Material architecture is also important 

- model for A ^ model for B (surface roughness and porosity are modeled) 
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. Outline 

Microscopic scale simulation of the ablation of fibrous materials 


1 . Models and simulation tool 

- Material models (A vs. B) 

- Studied Problem 

- Simulation tool : AMA 

2. Simulation and analysis 

- Simulated Problem 

- A vs. B: Moderate Thiele number 

- A vs. B: Small Thiele number 

3. Discussion 

- Effective Reactive Surface Area 

- Effective reactivity 

4. Conclusion 
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1 . Models 

Materials: 2 ideal low density carbon/phenolic ablators 



Preform: carbon fibers, random orientations 
Fibers: diameter (10 |jm), length (0.5 mm) 

Virgin mass fractions: carbon fiber (65%), Phenolic resin (35%) 
Overall density (Virgin : 280 kg/m 3 ; Pyrolyzed: 230 kg/m 3 ) 


Similarities 


Preform: carbon fibers, random orientations 
Fibers: diameter (10 pm), length (0.5 mm) 
Virgin mass fractions: carbon fiber (65%), Phenolic resin (35%) 
Overall density (Virgin : 280 kg/m 3 ; Pyrolyzed: 230 kg/m 3 ) 


Dense phenolic resin surrounding the fibers: 1 pm 

Virgin phenolic resin density: 1 200 kg/m 3 
Pyrolyzed phenolic resin density: 600 kg/m 3 
Overall porosity: 0.85 


Difference 


Low density pore-filling matrix 

Virgin phenolic resin density: 1 00 kg/m 3 
Pyrolyzed phenolic resin density: 50 kg/m 3 
Overall accessible porosity: 0; closed porosity 





1. Models and simulation tool 

Ablation model: Transport, Reaction, and Local Surface Recession 


• Problem studied: Isothermal oxidation of materials A and B in their charred form 


• Model 

Starting point : differential recession of a 
heterogeneous surface Sby gasification 


dS 


+ \.VS = 0 


dt / 

Local recession velocity 


V — —J Q n ; i= {fiber , matrix} 

/ 

Local ablation flux per surface unit 


J 



i = {fiber , matrix} 


Oxygen transport 


dX 

+ V • (-DVX ) + v • VX = 0 

dt s 


Local Diffusion Convection 

concentration (if any) 



k = Reactivity (m s -1 ) 
i? = Normal to the surface ( - ) 
v g = Pyrolysis gas velocity (m s _1 ) 
O = Solid molar volume (m 3 mol -1 ) 
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1. Models and simulation tool 

Flow regime in the pores of the material [1] : from Knudsen to continuum 


Just as an illustration, Knudsen number in the pores along the Stardust trajectory 



Stardust trajectory 
T-P conditions at TPS surface 
Time tram entry interface 
(1 arrow = t s) 

21 s : beginning of pyrolysis 
51 s : peak heating 
1-33 s : end of ablation 


Mean free path of 0 2 in air 

X = 3.2x10" 


- 5 T(inK) 


p {in Pa) 


Knudsen number 

/ 1 {mean free path) 
d {pore diameter) 


[1] J. Lachaud, I. Cozmuta, N. N. Mansour. Multiscale approach to ablation modeling 

of phenolic impregnated carbon ablators. Journal of Spacecraft and Rockets, accepted for publication. 
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1. Models and simulation tool 

Random Walk for reaction/diffusion & Triangle Marching Cube for surface recession 
Random Walk : Knudsen & Intermediate (classical) - Continuum regime (Brownian Motion) 


3D simulation tool : AMA^ 



Boundary 
condition 
(Buffer zone) 

Mass transfer 
(classical random 
waik - small steps 
following Maxwell- Boltzmann 
distribution - 
or Brownian motion 
method - large steps to lower 
computation lime in 
continuum regime) 

Ablation rate 
(Slicking event) 

Moving fluid/solrd 
interface 
(simplified 
Marching cube) 


J. Lachaud, G.L. Vignoles. A Brownian 
motion technique to simulate gasification 
and its application to C/C composite 
ablation. Computational Materials 
Science, 44(4), 2009, pp. 1 034-1041 . 
doi :1 0.101 6/j.commatsci.2008.07.01 5 
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1. Models and simulation tool 

Determination of the effective diffusion coefficient from first principle 


Illustration : path of a walker in a periodic cell (code AMA). 
The material in this illustration is anisotropic. 



The fibrous media is tortuous. This slows 
the diffusion process (collisions on the 
walls). We can determine the tortuosity 
factor as a function of the mean free path 
of the molecules; that is, as a function of 
the Knudsen number. 
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. Outline 

Microscopic scale simulation of the ablation of fibrous materials 


1 . Models and simulation tool 

- Material models (A vs. B) 

- Studied Problem 

- Simulation tool : AMA 

2. Simulation and analysis 

- Simulated Problem 

- A vs. B: Moderate Thiele number 

- A vs. B: Small Thiele number 

3. Discussion 

- Effective Reactive Surface Area 

- Effective reactivity 

4. Conclusion 
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2. Simulation and analysis 

Fiber preform : 1 D steady-state analysis (diffusion time « ablation time) 


• The concentration field is a function of the Thiele number 




X( coshlWL^ 


coshO 


Thiele number 
0 = L 


L s : material depth (m) 

D eff : effective diffusivity (m 2 /s) 
k f : fiber reactivity (m/s) 
s : specific surface (m 2 /m 3 ) 


^(z=o/^0“^ 

surface 


C(z)/C n 


Large Thiele number : ablation is a surface phenomenon, A = B 


Smaller Thiele number : ablation is a volume phenomenon, A ^ B ? 


bottom 



2. Simulation and analysis 

Oxidation of B at moderate Thiele Number, O = 40 


• Hypotheses 

- B is fully charred 

- No pyrolysis gas blowing 

- Air 

- O = 40, e. g. 

• T = 3360 K (isothermal) 

• P = 0.26 Atm 

• Pi = 32p m 

• oxidation by 0 2 
k m = 10 k f = 13.7 m/s 
[Drawin 1992, Lachaud 2007] 

N.B.: oxidation by O 

k f = 1 00 m/s [Park], O > 1000 


■ 


Oxidation by the top 



Periodic on 
the sides 




• Diffusion/reaction DSMC simulation 

- Physical time: 1 .2 s 

- Computational time : 2 days on 
a single core. 


12 


0.75 mrr\ 


2. Simulation and analysis 

Oxidation of B at moderate Thiele number (<t> = 40): 3 stages 


1) Matrix is removed 2) Fiber diameter 3) Overall surface 
first decreases recession 



t = 0.0 s 


0.4 s 


0.8 s 


1.2 s 


13 


ablation zone 



Avizcr 


2. Simulation and analysis 

Same conditions (<t> = 40): Comparison of A & B 



2. Simulation and analysis 

Same conditions (O = 40): Comparison of A & B 


Material A 


Material B 


Simulation time: 0 % 


5 


0.6 s 






Similar behaviors at high and moderate Thiele number 

i.e, high temperatures 
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2. Simulation and analysis 

Small Thiele number (O < 0.01 ); Material A 



t = 0.0 t.u. 


5 t.u. 


10 t.u. 


100 t.u. 


AtT=1000 K, 1 t.u. «4 minutes, oxidation by molecular oxygen. 
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2. Simulation and analysis 

Small Thiele number (O < 0.01 ): comparison of A and B 


Material B 


Simulation time*: 0 


5 


100 




Different behaviors at small Thiele numbers 
i.e., moderate temperatures 
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. Outline 

Microscopic scale simulation of the ablation of fibrous materials 


1 . Models and simulation tool 

- Material models (A vs. B) 

- Studied Problem 

- Simulation tool : AMA 

2. Simulation and analysis 

- Simulated Problem 

- A vs. B: Moderate Thiele number (0=40) 

- A vs. B: Small Thiele number (0=40) 

3. Discussion 

- Effective Reactive Surface Area 

- Effective reactivity model 

4. Conclusion 
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ERSA 


3. Discussion 

Effective Reactive Surface Area (ERSA). Illustration: fiber preform [1]. 


For a fiber preform, the effective reactivity is k eff = Versa kf(T) 


Versa 


ERSA 

GSA 


L 

J s f X(z)/X 0 dz 

z = 0 


T tanhd 
sL 
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+(i-s) 


lim ' s L + (1-f) 



l °Sl0 D eff 

Plotted for: k f =1 m/s; L=1cm; s=6x10 5 rrf/m 3 


L s : material depth (m) 

D eff : effective diffusivity (m 1 2 /s) 
k f : fiber reactivity (m/s) 
s : specific surface (m 2 /m 3 ) 

GSA: geometric surface area (m 2 ) 
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[1] J. Lachaud, I. Cozmuta, N. N. Mansour. Multiscale approach to ablation modeling 

of phenolic impregnated carbon ablators. Journal of Spacecraft and Rockets, accepted for publication. 
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3. Discussion 

Effective reactivity of material B at small Thiele number (O < 0.01) 



In this case, the effective 
reactivity can be estimated 
analytically under the 
following hypotheses: 

- No temperature gradient 

- No pyrolysis gas flux 

- The fiber diameter 
reduction is small 


We are currently working on the development of more elaborated models 
including 

- Temperature gradients, 

- Homogeneous finite-rate chemistry (ablation and pyrolysis gases) 


age) 
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m m j j 


) 


0.0007 


k 0 : reactivity at t=0 (m/s) 
s : specific surface (m 2 /m 3 ) 
h: matrix depth (m) 
e f : fiber volume fraction 
k f : fiber reactivity (m/s) 

£ m : matrix volume fraction 
k m : matrix reactivity (m/s) 
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. Conclusion and perspectives 


s New volume ablation model (oxidation only) 
o Large Thiele number : A* B 
(usually high temperatures) 
o Smaller Thiele number : B 

(usually moderate and low temperatures) 

s Finite-rate chemistry 

o The effective reactivity of a material is not 
only a function of the temperature. It is also a 
function of the ERSA; that is, of: 

- the geometric surface area available 
(depends on the ablation history); 

- the Thiele number; 

- the temperature gradient; 

- homogeneous reactions occurring with 
pyrolysis gases. 


s Perspective: Application to real re-entry conditions 
o Thermal gradients 

o Pyrolysis-ablation coupling (pyrolysis gases) 
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Appendix 



. Application to Stardust conditions 

Material B : Fit of post-flight density profile [1] 
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— Surface ablation 
(current models) 


(j) Stardust Post-Flight Analyses 
[Stack poo le, 2008] 


Volume ablation model 
(this study) 
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[1] J. Lachaud, I. Cozmuta, N. N. Mansour. Multiscale approach to ablation modeling 

of phenolic impregnated carbon ablators. Journal of Spacecraft and Rockets, accepted for publication. 
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. Current developments: chemistry 

Pyrolysis-ablation coupling. Main problem : Finite-rate chemistry. 


Shock 



CFD : 

-finite rate chemistry 
(homogeneous only) 

- diffusion 

- convection 


M aterial code : 

- ablation H-pyrolysis 
-^finite rate 
chemistry 
(homogeneous and 
heterogeneous) 

- pyrolysis : multi- 
laws 

- diffusion 

- convection 
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. Current developments: macroscopic scale model 

Comparison of A & B for the end of Stardust re-entry 



Surface ablation models 

Material B + Volume ablation model 

Material A + Volume ablation model 

1 D simulation of the evolution 
of the density profile of the 
char layer of A &B in the 
conditions encountered during 
the end of the reentry of 
Stardust 

Trajectory time 
90s < t < 130s 
Wall temperature 
1300 K > T > 850 K 

Temperatures provided by the 
TPS-ADP / Stardust post-flight 
analysis team. 

Hypothesis 
Thermal Conductivity 

^a=^b= kpicA 


z ( mm) 
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